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A search for muon neutrinos from Kaluza-Klein dark matter annihilations in the Sun has been 
performed with the 22-string configuration of the IceCube neutrino detector using data collected 
in 104.3 days of live-time in 2007. No excess over the expected atmospheric background has been 
observed. Upper limits have been obtained on the annihilation rate of captured lightest Kaluza- 
Klein particle (LKP) WIMPs in the Sun and converted to limits on the LKP-proton cross-sections 
for LKP masses in the range 250 - 3000 GeV. These results are the most stringent limits to date on 
LKP annihilation in the Sun. 



PACS numbers: 95.35.+d, 98.70.Sa, 96.50.S-, 96.50.Vg 

In a recent work [l[ , we presented the result of a search 
for neutralino dark matter accumulated in the center of 
the Sun with the 22-string configuration of the IceCube 
detector. In this letter we extend the search to an alter- 
native dark matter candidate, Kaluza-Klein (KK) parti- 
cles, arising from theories with extra spacetime dimen- 
sions. In the simplest framework of universal extra di- 
mensions (UED) 0], there is a single extra dimension of 
size R ~ O (TeV^ 1 ) compactified on an S 1 jZi orb- 
ifold. Within minimal UED theories, the first excitation 
of the hypercharge gauge boson, is generally the 

lightest KK particle (LKP). It is often denoted as the 
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KK 'photon', 7W, since the effective first KK-level Wein- 
berg angle of the mass matrix is very small, and therefore 
£?W can also be described as a mass eigenstate p]- KK- 
parity conservation, affiliated with extra-dimensional mo- 
mentum conservation, leads to the stability of the LKP, 
which makes it a viable DM candidate. There are also 
other possible natural choices for LKP candidates within 
UED, like the KK 'graviton', the KK 'neutrino' or the 
Z^-boson that may constitute viable DM candidates. 
They are not considered here. Instead, we focus on the 
most promising KKDM prospect in terms of indirect de- 
tection expectations, the KK 'photon'. 
Accelerator measurements constrain the lower bound for 
the mass of the LKP, m 7 (i), at 300 GeV [3j. The upper 
bound is limited to a few TeV in order to not exceed the 
observed DM relic density and overdose the Universe. 
We here consider UED models with five spacetime di- 
mensions characterized by two parameters: the LKP 
mass, m 7 (i>, and the mass splitting A g (i) = (m ? (i> — 
m 7 (i))/m 7 (i), where m q (i) is the mass of the first KK 
quark excitation, as discussed in @, |3> IE • 
As a possible dark matter component of the halo, LKPs 
can become gravitationally trapped in massive celestial 
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TABLE I: LKP annihilation branching ratios for two values 
of A 3 (i) [7j. Ratios are not summed over generations. Chan- 
nels within parenthesis give negligible contribution to a neu- 
trino flux from the Sun. The Higgs-field annihilation channel, 
marked with ' , is neglected, due to large uncertainty and small 
contribution to the neutrino flux. 
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FIG. 1: Comparison of simulated muon spectra from LKP, 
7^ 1 - ) , and neutralino, Xii annihilations observed in IceCube, 
for two WIMP masses, 250 and 3000 GeV, representing the 
boundries of the investigated LKP model space. 



bodies like the Sun, accumulating to high DM densities 
that can exceed the mean galactic density by several or- 
ders of magnitude in the object's core. Since the LKP 
is a boson, pair-wise annihilation is dominated by s-wave 
processes, creating standard model particles whose de- 
cay chains produce neutrinos in the GeV - TeV range. 
The branching ratios for the LKP annihilation channels 
of interest are given in Table Q] for two values of A (i) 
7}. The neutrinos may escape the Sun and reach Earth. 
The search presented here aims at detecting LKP anni- 
hilations indirectly by observing an excess of such high 
energy neutrinos from the direction of the Sun. Despite 
the existence of various limits on neutralino induced neu- 
trino fluxes from the Sun [l|, la La llIA > no correspond- 
ing limits for LKP annihilations have been previously 
reported. 

For the results presented here, we use the same data 
set, 104.3 days livetime taken with the 22-string configu- 
ration of IceCube in 2007, and the same analysis cuts as 
presented in [l[ . This is justified since the signature of the 
expected signal at the detector is very similar for the LKP 
and neutralinos, considering the hardest x5~ annmna tion 
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FIG. 2: Lines showing the effective area (left scale) for the fi- 
nal event selection as function of neutrino energy in the range 
50-1000 GeV, for muon neutrinos (solid line) and antineutri- 
nos (dashed line) from the direction of the Sun. The result is 
an average over the austral winter. Systematic effects are in- 
cluded at the la level, and statistical uncertainty of the same 
level are shown with error bars. Also shown is the median 
angular error with la error bars (squares, right scale). 



channel into W + W~ . The neutrino spectrum from an- 
nihilations of a LKP of a given mass in the center of the 
Sun is considerably harder than that of a neutralino of 
the same mass. However, oscillations and energy losses 
of the neutrinos on their way out of the Sun, like neutral 
current (NC) scattering, absorption and ^-regeneration, 
smear out the energy spectra in a way that makes them 
comparable at Earth. Figure [1] shows an example of how 
the resulting muon spectra at the detector compare for a 
selected choice of neutralino and LKP masses at 250 and 
3000 GeV. The analysis strategy used in [l[ is therefore 
already optimized for the search of KK dark matter. 

We simulated LKP annihilations in the Sun 
using WimpSim [12] for LKP masses 



m 7 (i) 



250,500,700,900,1100,1500,3000 GeV. We used 
A q (i) = with annihilation branching ratios from 
Table [IJ Since A ? (i) > results in an increased neutrino 
flux due to the importance of the contributions from the 
t + t~ and the direct neutrino channels, the choice of 
A„(i) = leads to a conservative limit. The background 
in the search for neutrinos from the Sun comes from air 
showers created in cosmic ray interactions in the atmo- 
sphere. The showers cause downward going atmospheric 
muon events, triggering the detector at ~ 500 Hz, and 
atmospheric muon neutrino events, triggering at ~ 4 
mHz. When the Sun is below the horizon, the neutrino 
signal can be distinguished from the atmospheric muon 
background by selecting events with upward-going 
reconstructed tracks. 

Atmospheric muon and neutrino background events 
were also generated [l3|, [lij . The pr opagat ion of muons 
and photons in the ice was simulated [la . Il6l ] taking mea- 
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TABLE II: Upper limits on the number of signal events /j, a , the conversion rate r„_ M , the LKP annihilation rate in the Sun Ta, 
the muon flux "I^, and the LKP-proton scattering cross-sections (spin-independent, <r SI , and spin-dependent, <r SD ), at the 90% 
confidence level including systematic errors. The sensitivity "!> M (see text) is shown for comparison. Also shown is the median 
angular error 0, the mean muon energy (E M ), the effective volume V e g, and the effective area A e g. 
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sured ice properties into account [l7j |. 

The events had to pass several selection criteria as de- 
scribed in [l[ in order to reduce the content of atmo- 
spheric muon events. As a compromise between signal 
efficiency and background rejection, it was required that 
more than half of the events in the final data sample 
were neutrino-induced. The observables used describe 
the quality of the track reconstructions and the geome- 
try and time evolution of the hit pattern in the detector, 
and they were required to be well reproduced in simu- 
lations. The event selection consisted first in a series of 
unidimensional cuts on the selected event variables, and a 
final step that used two Support Vector Machines (SVM). 
The SVMs were trained with simulated signal, and a set 
of experimental data, recorded in December 2007 and not 
used in this analysis since the Sun was above the horizon, 
was taken as background. A final sample was then de- 
fined from a cut on the combined two SVM output values, 
Qi x Qi (see figure 1 in [l[). The analysis was performed 
in a blind manner such that the azimuth of the Sun is 
unknown until the selection criteria were finalized. 

The systematic uncertainties on the effective volume, 
V e g, defined as the equivalent detector volume with 100% 
selection efficiency, are the same as the ones calculated 
in the WIMP analysis in [l[, and are dominated by the 
uncertainties in photon propagation in the ice and the 
absolute DOM efficiency. They range from ±19% for the 
highest m 7 (i> to ±26% for the lowest m 7 <i) [la]. From 
the final event selection of the signal simulation we addi- 
tionally derive the effective area for muon neutrinos from 
the direction of the Sun as a function of neutrino energy, 
see Fig. [2] Also shown in the figure is the median angular 
error, the median of the angle between the reconstructed 
muon and the neutrino direction, O. The result includes 
systematic uncertainties and is an average over the aus- 
tral winter, during which the Sun is below the horizon. 

For the LKP signal models we then calculated the effec- 
tive volume and, based on the distribution of the recon- 
structed angle to the Sun 'J, we constructed confidence 
intervals at the 90% confidence level using the method 
outlined in [l|: to evaluate the signal content in the fi- 
nal event sample, hypothesis testing was done based on 
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FIG. 3: Cosine of the angle between the reconstructed track 
and the direction of the Sun, for data (squares) with one 
standard deviation error bars, and the atmospheric back- 
ground expectation from atmospheric muons and neutrinos 
(dashed line). Also shown is a simulated signal (m (i) = 900 
GeV) scaled to fi s — 7.0 events (see Table ITT)| . 



"J, the angle between the reconstructed track and the 
direction of the Sun. From simulations we find f s (*&), 
the probability distribution of W for the signal. By ran- 
domizing the azimuth angle in the final event sample of 
experimental data, /b(4'), the equivalent probability dis- 
tribution is found for background. Defining £ = 
from the number of signal events fi s and the observed 
number of events n bs, we form the combined probabil- 
ity density = £ • /„(¥) + (1 - £) • / b (tf). Based 

on £bestj the non-negative signal content that maximizes 
the likelihood, we form the logarithm of the likelihood 



ratio R(0 = log(n-=?° bs f^jl) ) UM- Comparing this 
with a -Rtcst(C) distribution of a large number of pseudo- 
experiments with n bs events taken from /f (^), we con- 
struct the confidence interval on £ at significance a as 
R(£um) = RfcstiZhm), where P(R tcs t > R? cst ) = I - a. 
No excess of events from the Sun above the background 
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FIG. 4: Limits on the muon flux from LKP annihilations 
in the Sun including systematic errors (squares), compared 



to the theoretically allowed region of m (i) and A 



The 



regions corresponding to A (ij = 0.01 and A (i) =0.1 are 
marked with black lines. The region below m^(i) = 300 GeV 
is excluded by collider experiments 



expectation was found in the search (£t>cst = 0) . The ob- 
served number of events as a function of the angle to the 
Sun, 'J, is compared to the atmospheric background ex- 
pectation in Fig. [3] From the upper limits on the number 
of signal events fi s we calculate the limit on the neutrino 
to muon conversion rate r„_>,, = ./*' . , for the livetime 
t. Using the signal simulation [12j |. we can convert this 
rate to a limit on the LKP annihilation rate in the Sun, 
Ta, see Table ILT1 Results from different experiments are 
commonly compared by calculating the limit on the muon 
flux above 1 GeV, $ M , which is also given in Table |TT] to- 
gether with the sensitivity, $ M , the median limit obtained 
from simulations with no signal. 

The flux limit is shown in Fig. 0] together with the the- 
oretically allowed flux region, derived from Refs. [f| 0] 
with the use of DarkSUSY |20| . We have here approxi- 
mated the branching ratios for the regions of A g (i) = 0.01 
and A ? (i) = 0.1 with those of A„(i) = and A 9 <i> = 0.14, 
respectively, as given in Table ffl The limits on the an- 
nihilation rate can be converted into limits on the spin- 
dependent, er SD , and spin- independent, er SI , LKP-proton 
cross-sections, allowing a comparison with the results of 
direct search experiments. Since capture in the Sun is 
dominated by cr SD , indirect searches are expected to be 
competitive in setting limits on this quantity. Assum- 
ing equilibrium between the capture and annihilation 
rates in the Sun, the annihilation rate is directly pro- 
portional to the cross-section. A conservative limit on 
cr SD is found by setting a SI to zero, and vice versa. We 
have used the method described in Ref. [2l[ to perform 
the conversion. The results are shown in Table |TTJ We 
assumed a local WIMP density of 0.3 GeV/cm 3 , and a 
Maxwellian WIMP velocity distribution with a disper- 



FIG. 5: Limits on the LKP-proton SD scattering cross-section 
(squares) adjusted for systematic effects, compared with lim- 
its from direct detection experiments [22 . 123 . 24] . Theoreti- 
cally predicted cross-sections are indicated by the green area 
5]. The regions corresponding to A ? (i) = 0.01,0.1,0.5 are 
marked with black lines. The region below m^(i) = 300 GeV 
is excluded by collider experiments and the upper bound 
on m 7 (i) corresponds to the overclosure limit for each individ- 
ual LKP model [J. The lighter blue region is allowed when 
considering 0.05 < !2odm'i 2 < 0.20, and the darker blue re- 
gion corresponds to the preferred la WMAP (5 year) relic 
density 0.1037 < Ucdm/i 2 < 0.1161 H|. 



sion of 270 km/s. Planetary effects on the capture were 
neglected (2(| . Figure [5] shows the limits on er SD , as ob- 
tained with the 22-string configuration of IceCube com- 
pared with other bounds [22|, HH [13] , and the KK model 
space. The theoretical model space (green area) is plot- 
ted for different predictions for the mass splitting A ? <i) . 
The blue regions indicate the overlap regions with two 
different Ocdm intervals, whereas the narrow dark blue 
region corresponds to the preferred WMAP lc-region for 
CDM. The upper bound on m^i) , derived from the over- 
closure limit for each individual LKP model [4|, varies 
with different values of A ? (i) and increases remarkably 
for models with A 9 (i) < 0.1. This is due to additional 
coannihilation effects, arising for degenerate LKP mod- 
els @. 

In conclusion, we have presented the first limits on 
LKP annihilations in the Sun. We also derived the 
most stringent limits on the spin-dependent LKP-proton 
cross sections in the non-excluded LKP mass regions 
(300GeV < m 7 (i> < 3TcV), improving existing limits 
by more than two orders of magnitude and excluding 
some viable LKP models. The full IceCube detec- 



tor with the DeepCore extension [27j is expected to 
test most LKP models within the allowed region for 
0.05 < Ocdm/j 2 < 0.20, shown in Fig. [5] 
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